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Abstract 

The superfluid density is a fundamental quantity describing the response to a rotation as well as 
in two-fluid collisional hydrodynamics. We present extensive calculations of the superfluid density 
p s in the BCS-BEC crossover regime of a uniform superfluid Fermi gas at finite temperatures. We 
include strong-coupling or fluctuation effects on these quantities within a Gaussian approxima- 
tion. We also incorporate the same fluctuation effects into the BCS single-particle excitations 
described by the superfluid order parameter A and Fermi chemical potential u, using the Nozieres 
and Schmitt-Rink (NSR) approximation. This treatment is shown to be necessary for consistent 
treatment of p s over the entire BCS-BEC crossover. We also calculate the condensate fraction iV c 
as a function of the temperature, a quantity which is quite different from the superfluid density 
p s . We show that the mean- field expression for the condensate fraction iV c is a good approxima- 
tion even in the strong-coupling BEC regime. Our numerical results show how p s and N c depend 
on temperature, from the weak-coupling BCS region to the BEC region of tightly-bound Cooper 



pair molecules. In a companion paper by the authors (cond-mat/0609187), we derive an equiv- 
alent expression for p s from the thermodynamic potential, which exhibits the role of the pairing 
fluctuations in a more explicit manner. 

PACS numbers: 03.75. Ss, 03.75.Kk, 03.70.+k 
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I. INTRODUCTION 



In the last few years, the BCS-BEC crossover in two-component Fermi superfluids has 
become a central topic in ultracold atom physics Q, 0, 3, I^J. This crossover is of especial 
interest since the superfluidity continuously changes from the weak-coupling BCS-type to 
the Bose-Einstein condensation (BEC) of tightly bound Cooper pairs, as one increases the 
strength of a pairing interaction [5| . Thus the BCS-BEC crossover enables us to study fermion 
superfluidity and boson superfluidity in a unified manner. 

The superfluid density p s is a fundamental quantity which describes the response of 
a superfluid which arises from a BEc((|. The superfluid density was first introduced by 
Landau as part of the two-fluid theory of superfluid 4 HeJ3|. At T = 0, the value of p s always 
equals the total carrier density n. (In the BCS-BEC crossover, n is the total number density 
of fermions.) This property is satisfied in both the Fermi superfluids and Bose superfluids, 
irrespective of the strength of the interaction between particles. This is quite different from 
what is called the condensate fraction N c , which describes the number of the Bose-condensed 



particlesll , [8|] . For example, in superfluid 4 He, only about 10% of atoms are Bose-condensed 
even at T = 0, due to the strong repulsion between the 4 He atoms [For a review, see Ref. 

n 

|9j.]. In contrast, all the atoms contribute to the superfluid density at T = 0, namely 
Ps (T = 0) = n. 

In a companion paper, we have discussed some analytical results for the superfluid density 
p s in the BCS-BEC crossover regime of a uniform superfluid Fermi gasQ|. Going past the 
weak-coupling BCS theory, we derived an expression for p s in the Gaussian fluctuation level 
in terms of the fluctuations in the Cooper channel. The resulting expression for the normal 
fluid density p n = n — p s is given by the sum of the usual BCS normal fluid density p^ and a 
bosonic fluctuation contribution p^. While the superfluid density from fermions dominates in 
the weak-coupling BCS regime, the bosonic fluctuation contribution p 1 ^ becomes dominant in 
the strong-coupling BEC regime. Since p^ is absent in the mean-field BCS theory, inclusion 
of fluctuations in the Cooper channel is clearly essential in considering the superfluid density 
in the BCS-BEC crossover. 

In Ref. [10|, our expression for p s was obtained using the thermodynamic potential in 
the presence of a superfluid flow. In the present paper, we derive a second expression for 
p s by calculating the effect of pairing fluctuations on the single-particle Green's function 
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with a supercurrent. In this paper, we use this expression to numerically calculate p s in 
the entire BCS-BEC crossover regime at finite temperatures. However, it can be shown 
ihat this expression for p s (given in Sec. Ill) is equivalent to the result derived in Ref. 
lp| . In calculating the superfluid order parameter A as well as Fermi chemical potential 
fi, we include the effect of the pairing fluctuations following the approach given in Ref. J4J]. 
This self-consistent treatment of A and \x is crucial in calculating function of the 

temperature in the BCS-BEC crossover. 

Besides the superfluid density, we also calculate the condensate fraction N c describing the 
number of Bose-condensed particlesjjj. N c is of special interest in superfluid Fermi gases, 
since it can be observed experimentally. Indeed, a finite value of N c is the signature of the 
BCS-BEC superfluid phase j5|. In this paper, we show that strong-coupling pair fluctuations 
have little effect on N c in the BCS-BEC crossover. We note that N c has been recently 
calculated at T = within a simple mean-field BCS approach [ill], and using Monte Carlo 
techniques}"!^!]. In this paper, we present detailed results for N c at finite temperatures in the 
BCS-BEC crossover, based on the NSR theory of fluctuations |2|. 

The present paper is organized as follows. In Sec. II, the BCS Green's functions are solved 
numerically for the superfluid order parameter A and chemical potential fi self-consistently. 
This is done for the entire BCS-BEC crossover region and at finite temperatures. In Sec. Ill, 
we calculate the superfluid density treating the strong-coupling pair fluctuation effects 
within a Gaussian approximation 2, SO, 13- Numerical results for p s are presented in Sec. 
IV. In Sec. V, we define and calculate the condensate fraction N c . 

Throughout this paper, we take Ti = ks = 1. We also set the volume V — 1, so that the 
number of atoms N and the number density n are the same. 



II. BCS-BEC CROSSOVER IN THE SUPERFLUID PHASE 

In superfluid Fermi gases, current experiments make use of a broad Feshbach resonance 
to tune the magnitude and sign of the pairing interaction^. In this case, the superfluid 
properties can be studied by using the single-channel BCS model, described by the Hamil- 
tonian, 

H = Y,£p C U C P° - U E C p+qT C i'-q| C P'! C pT- (2- 1 ) 

P' 7 p,p',q 
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Here, c pcr is a creation operator of a Fermi atom with pseudo-spin a =],[ (which describe the 

two atomic hyperfine states.). The Fermi atoms have kinetic energy £ p = e p — \i = p 2 /2m—{i, 

measured from the Fermi chemical potential /i. —U describes a pairing interaction between 

different Fermi atoms. The magnitude and sign of U can be tuned using the Feshbach 

resonance using an external magnetic field. The weak-coupling BCS limit corresponds to 

U — > +0. We only consider a uniform gas in this paper. 

Nozieres and Schmitt-Rink (NSR) first discussed the BCS-BEC crossover behavior of 

Q 

the Hamiltonian in Eq. (|2.1|) to determine T c |2j, which we brieflv summarize (see also Ref. 
jl^). based on a Gaussian approximation for pair fluctuations sj]. The NSR theory gives 
two coupled equations, one for T c , and one giving \i as a function of N. The equation for T c , 
given by the Thouless criterion, has the same form as the weak-coupling BCS gap equation 

at T c , 

l = ^E^tanhfe P , (2.2) 

where (3 — 1/T. In the weak-coupling BCS theory, one finds that the equation of state gives 
\x = £p, where is the Fermi energy. However, this result for // is not valid in the BCS-BEC 
crossover regime where we must include strong-coupling effects due to pairing fluctuations. 
In the NSR theory, this deviation is determined by solving the equation for the number of 



fermions 



15|, within a t- matrix approximation in the particle-particle channel, 



n = 2J2 /(£ P ) - ~ E M 1 - ^n(q> iv n )\. (2.3) 

Here /(£ p ) is the Fermi distribution function, and n(q, iv n ) describes non-interacting fluc- 
tuations in the Cooper channel, 

n(q, iu n ) = E 1 : /(ep t L q/ : ) " /(ep : q/2) - (2-4) 

p ^p+q/2 + ?p+q/2 ~ 

Here, v n is the boson Matsubara frequency describing bosonic fluctuations. The second term 
in Eq. (|2.3J) gives the number of Cooper pairs which are formed, and which become stable 
as we go through the Feshbach resonance. In the extreme BEC limit, all the fermions form 
such stable molecules. 



The NSR theory to the superfluid phase has been studied below T c j4j, |13|]. We briefly 
review this approach. It is convenient to introduce the two-component Nambu field 
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operator [ly, III I 



(2.5) 



18], 



In this standard Nambu representation, Eq. (|2.1|) can be rewritten as 

# = ^ + E^p + E - An]* p - j + P2, q p 2 ,- q ]. (2.6) 

17 p p q 

The mean-field BCS approximation is described by the first three terms of Eq. (|2.6|) . with 
the 2 x 2-matrix Nambu-Gor'kov propagator [if]] 

1 



G (p,iu n 



+ £p7~3 - ATI 



llUJ r 



\2 



(2.7) 



p 

where u; m is the fermion Matsubara frequency associated with the fermions and r,- are the 



Pauli operators. E p = «/£ 2 + A 2 describes Bogoliubov single-particle excitations, associated 
with the breaking of a Cooper pair. The order parameter A = U Z) P ( c P t c - P |) (taken to be 
real) corresponds to an off-diagonal mean-field and is related to the thermal Green's function 
G (p,iw m ) as 

A = ^ E GffW m ). (2.8) 
Substituting Eq. ()2.7|) into Eq. (|2.8|1 . we obtain the well-known gap equation, 

A = ^E^tanh^. (2.9) 

Equation ()2.9j) is an implicit equation determining A as a function of fi and T. In the 

path integral formalism, the result in Eq. (12.91) emerges as the saddle point solution of Eq. 
rn 

(I2~K1)[3]. Equation (Q reduces to Eq. when A = 0, which defines 
Pj,i U = 1> 2) represent the generalized density operators, defined by(l9j] 

p 

Writing out the operators, one sees that pi ;q and p 2 , q describe amplitude and phase fluctua- 
tions of the Cooper-pair order parameter A, respectively. As shown in Eq. ([2.6)1 . the inter- 
action separates into amplitude fluctuations (pi,qPi,-q) and phase fluctuations (/?2,q/°2,-q)- 
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FIG. 1: Fluctuation contribution to the thermodynamic potential £1. The solid line describes the 
single particle matrix Green's function Go, and the dashed line represents the pairing interaction 
—U. II?. = 1,2) are the generalized correlation functions, given by Eqs. (|A5|) - (|A7|) . 

The chemical potential \x is determined by the equation for the total number of particles 
N. This equation is conveniently obtained from the thermodynamic identity N = —dtt/dfi. 
We calculate the thermodynamic potential Q perturbatively, taking into account the in- 
teraction in Eq. (12.61) . In the Nambu representation, the fluctuation contribution to Q is 

h nri 

diagrammatically described by Fig. [T] As discussed in Refs. |4j,|l3|, one finds [See Appendix 
A.] 

N = N° F - -LA J2 lndet[l + E7S(q,zV n )], (2.11) 

where the first term 

^ = E[l-#tanh§E p ] (2.12) 

comes from the mean-field thermodynamic potential Q F describing BCS Fermi quasiparti- 
cles. The second term in Eq. (|2.11|) describes the thermodynamic potential contribution Q B 
from bosonic collective pair fluctuations. H is a 2 x 2-matrix correlation function denned in 
Eq. ()A3)1 . We note that the //-derivative in Eq. ()2.11|) only acts on the chemical potential in 
the kinetic energy £ p . Equation (|2.11|) reduces to the NSR result in Eq. (|2.3|) at T c [where 
A = and Hn(q, iis n ) = H 2 2(q, —iv n ) — ~ n(q, iv n )]- See Appendix A for more details. 

We solve Eq. ()2.9|) together with Eq. (|2.11|) to give A and \i self-consistently. As usual, 
we need to introduce a high energy cutoff u c in these coupled equations. This cutoff can be 
formally eliminated by introducing the two-body s-wave scattering length a s ^, 

^ EE ^_ (2.13) 

Using a s in place of U, one can rewrite Eqs. (|2.9|) and (j2.11j) in the form 

47ra, x ^r 1 , 8 1 i s 

1 = -^f^2E^2 E »-2r} (2 - 14) 





FIG. 2: (a) Off-diagonal mean-field A, and (b) Fermi chemical potential \x in the BCS-BEC 
crossover. The pairing interaction is measured in terms of the inverse of the two-body scattering 
length a s , normalized by the Fermi momentum fcp. In these panels, the dotted line shows T c as a 
function of (kpag) -1 . In the strong-coupling regime, the apparent first order behavior of the phase 
transition is an artifact of the NSR Gaussian treatment of pairing fluctuations (see text). 

" = A? - ~ E " ^(4*0 + (2-15) 

where the momentum sums are now no longer divergent. The weak-coupling BCS regime and 
the strong-coupling BEC regime are, respectively, given by (/cf^s) -1 ^ — 1 and (k^a^ 1 > 1 
(where kp is the Fermi momentum). The region —1 < (kpa s ) _1 < 1 is referred to as the 
"crossover regime." In Appendix B, we show that Eqs. (|2.14J1 and (12.15)1 are the limiting 
results obtained from a two-channel coupled fermion-boson model|4] when the Feshbach 
resonance is broad. 

Figure 121 shows our self-consistent solutions of the coupled equations (|2.14|) and (|2.15|) in 
the BCS-BEC crossover at finite temperatures. These reproduce the NSR results for T^, 
as shown explicitly in Fig. El 

When one enters the crossover region, Fig. 0] shows that the order parameter A deviates 
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FIG. 3: (a) Superfluid phase transition temperature T c , and (b) chemical potential /i(T = T c ) 
in the BCS-BEC crossover. In panel (b), fi at T = is also shown. The curve 'BEC gives the 
strong-coupling BEC limit, where one finds \i = —l/2ma 2 s . 
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FIG. 4: Calculated values of the superfluid order parameter A as a function of temperature. 'BCS' 
labels the weak-coupling BCS limit. The bendover near T c is an artifact of our NSR Gaussian 
treatment of fluctuations. 
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FIG. 5: Fermi chemical potential u as a function of temperature. For comparison, in the unitarity 



limit [too 



theory 



0], MC results 
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give ju/ep = 0.44 ~ 0.49, and an improved version of NSR 



2lJ gives n/sp = 0.4 ~ 0.47 just below T c 



from the weak-coupling BCS result. In the strong-coupling BEC regime, although the super- 
fluid phase transition approaches the value T c = 0.218Tpj3|, A(T = 0) continues to increase. 
As a result, the ratio 2A(T = 0)/T c in the BEC regime is larger than the weak-coupling BCS 
universal constant 2A(T = 0)/T c = 3.54. We recall that on the BEC side of the crossover, 
where /x is negative, the energy gap is not equal to A|4J. I13I|. 

The chemical potential is strongly affected by fluctuations in the Cooper channel, and 
becomes negative in the BEC regime, as shown in Figs. Efb) and dfb) . In the strong- 
coupling BEC regime, /x approaches /x = — l/2ma? j20l| . Although /x strongly depends on the 
magnitude of the interaction, Fig. El shows that the temperature dependence of /x is very 
weak in the entire BCS-BEC crossover. Our results in the unitarity limit are in quite good 
agreement with quantum Monte Carlo simulation jl^. as well as a more self-consistent 



version of NSR 



2l|. 



In Fig. El the apparent first-order phase transition in the BEC regime is an artifact of 
the approximate NSR theory we are using. The reason is as follows. In the BEC regime 
(where /x — the single-particle BCS excitations E p ■ 



/£p + A 2 have a large energy 
gap given by E g = \/ fi 2 + A 2 ~ |/x|. This energy gap still exists at T c , where A = 0. In this 
regime, we can set tanh(/3£ p /2) ~ 1 in Eq. fl23U) and f(E p ) ~ in Eqs. (fA5]) -(fA7 |) . Then, 
Eq. ()2.14j) reduces to the expression /x = — l/2ma%, and Eq. (|2.15jl becomes, by expanding 
the correlation functions Hy(q, iu n ) in powers of q and %v n [For the details, see a similar 
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calculation given in Ref. 

N 1 



iSu n 



= N c0 + - 2J — coth -w q - 1 

l q o» q 2 

= iV c0 + iV d . (2.16) 
Here, £_* = g 2 /2M ( M = 2m ) and 



u*=yJeB(eB + 2U M Na>) (2.17) 

is the Bogoliubov excitation spectrum in an interacting gas of Bose molecules. D in Eq. 
(I2.16|) is the Bose Green's function, describing Bogoliubov excitations 

(ivjl-u,* ' (2 - 18) 

where Um = ^ttcim/M (where in our theory and M = 2m) is the effective s-wave 

repulsive interaction between Cooper pairs. iV c o is given by 

A 2 \/2m 3 / 2 A 2 

N c o = J2jH= A • (2-19) 

p 4 Sp 167T W |yU| 



In Sec. IV, we prove that A^ c0 as defined in Eq. ()2.19|) corresponds precisely with the formal 
definition for the condensate fraction in a Fermi superfluid in the BEC limit. Na as defined 
in Eq. (J2.16|) is the number of molecules which are not Bose- condensed, again in the BEC 
limit. 

Equations (|2.16J) and (J2.17|) show that the strong-coupling BEC limit corresponds to the 
Popov approximation for a weakly-interacting molecular Bose gas [24]. As is well known (see, 



• 3; 



the Popov approximation gives a spurious first-order phase transition 



for example, Ref. 

at T c . This is the origin of the bendover or first-order phase transition evident in Fig. |2[26] 
and other figures in this paper. It is well known how to overcome this problem, na mely 
one has to include many-body renormalization effects due to the interaction Um 2Jj- l27 |. 



Including such higher order corrections past the NSR Gaussian fluctuations considered in 
this paper is also crucial for determining the correct value of effective interaction Um- The 
molecular scattering length clm = 2a s which is obtained in Eq. (|2.16[) is characteristic of the 
NSR treatment of fluctuations The correct result cim = 0.6a s 28l requires going past the 
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Gaussian approximation 



27 



291 ]. In this paper, in contrast, we only treat pairing fluctuations 



within NSR. However, within this approximation, we calculate the superfluid density (and 
condensate fraction N c in Sec V) in a consistent manner. 

As discussed in Ref. (25 1, the Popov approximation becomes invalid in the small region 



close to T c given by 



5t = ^— — < (-^) 1/3 (k F a M ) = 0.2Q(k F a M ). (2.20) 



Although we plot numerical results in the present paper in the whole temperature region for 
completeness, we emphasize that the restriction in Eq. ()2.20j) also holds in the BEC regime. 
We note that the region defined in Eq. (|2.20|) becomes narrow as one enters deeper into the 
BEC regime, simply because the molecular scattering length oc a s becomes small. Thus, 
one obtains 5t < 0.26 at (/cfo s ) _1 = 2 (the case shown in Fig. 01 for example) but 5t < 0.1 
at (/cpQs) 1 — 5. As Fig. 0] shows, the bendover occurs over an increasingly small region as 
we go deeper in the BEC region. 

Although Eq. ()2.16j) was obtained in the strong-coupling BEC regime, we note that 
the condensate fraction N c in Eq. (j2.19j) is the mean-field approximation N c for a Fermi 
superfluid. This is consistent with the result found in Sec. V that the mean-field expression 
for the condensate fraction N c is a good approximation even in the strong-coupling BEC 
regime (at least within NSR). The number of molecules in the non- condensate given in 
Eq. (|2.16|) . in contrast, is due to the pairing fluctuations. This is discussed in more detail 
in Sec. V. 



III. SUPERFLUID DENSITY AND THE SINGLE-PARTICLE GREEN'S FUNC- 
TION 

In Ref. JjLo|], our discussion of the superfluid density p s started from the thermodynamic 
potential Q(v s ) in the presence of an imposed superfluid velocity (or phase twist) v s . This 
is given by a simple generalization of Q given in Eq. (jA2J) . which is for the case v s = 0. 
Here we give an alternative formulation of p s in terms of how the single-particle Green's 
function is altered in the presence of a supercurrent. Our numerical calculations in Sec. IV 
are based on this expression, but it can be proven to be equivalent to the one discussed 
in Ref. [10]. The expression for p s discussed in Ref. jl0| is convenient for understanding 

11 



the role of collective pairing fluctuations. The result we obtain in this section gives further 
insight and is convenient for numerical calculations. 

When a supercurrent flows in the z-direction with the superfluid velocity v s = Q z /2m, 
the supercurrent density J z is given by 

J z = Y.— ( c P aCp, CT ) = nv s + - Yl — Tr ^(P>^m)]- (3- 1 ) 

P,(T m P p,UJm m 

Here g(p, iuj m ) is the 2 x 2- matrix single-particle thermal Green's function in the presence of 
v s . When the second term in Eq. (|3.1|) is expanded to 0(v s ), it can be written as J z = p s v s , 
where the superfluid density p s is defined by 

2 d 

Ps = U+ 3 ^ Pz dd~ Tr ^ ( " P ' iuJm ^Q*~*° = n ~ Pn- (3-2) 

The second line defines the normal fluid density p n . 

We first evaluate the second term in Eq. ()3.2j) in the lowest order mean-field approxima- 
tion. In this case, the supercurrent state is described by the BCS-type Hamiltonian 

~ 2 

H = Y,J cfrV»J(r)[^-^(r)- / dr[A(v)^\(r)4(r) + h.c.}, (3.3) 

where ?/v(r) is a fermion field operator, and 

A(r) = Ae jQ r , Q = (0, 0, Q z ) (3.4) 

is the superfluid order parameter for the current carrying state. In momentum space, Eq. 
(|3.3|) can be written as 



H = E^P C P,a C P,-- A E[ C P+QAT C -p+Q/2,| + 
Per p 

= E^U^s + ^p-Anl^p, (3.5) 



p 



where constant terms have been left out [see Eq. (J2.6|) ]. The effects of the supercurrent v s 
appear in the Doppler shift term, a p = Q • p/2m, and in £ p = e p — p, with p = p — Q 2 /8m. 
However, we do not have to take this dependence of £ p on p into account in calculating p n 
in Eq. ()3.2|) because it is of order of 0(vf). \l/ p is the two-component Nambu field in the 
supercurrent state, 

/ 



^p 



c p+Q/2,T 
V C - P +Q/2,| 
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(3.6) 



The BCS-Gor'kov Hamiltonian in Eq. ()3.5|) is described by the 2 x 2-matrix single-particle 
thermal Green's function with v s 7^ which has the well-known form [compare with Eq. 



g (p,ico m ) = - fdre 1 ^ (T r {* p (r)* P (0)}> = * — . (3.7) 

Jo (iuj m - otpj - 4 p t 3 + An 



Substituting Eq. (|3.7|) into Eq. (|3.2|) . one obtains the mean- field result 

^-ivW- (3 - 8) 

This is the expected Landau expression for the normal fluid density due to Fermi quasipar- 
ticles with the BCS spectrum E p when v s = 0. This expression gives the Fermi contribution 

to the total normal fluid density p n throughout the BCS-BEC crossover, and is not 
limited to the weak-coupling BCS limit. See, Ref. for additional discussion. 

In addition to p^ as given by Eq. ([3.8)1 . the correction to g to first order in v s gives 
rise to an additional fluctuation contribution p 1 ^ to the normal fluid density. This should be 
consistent with the number equation in Eq. (j2.11|) . for the v s = state, 

N = + 4 £ Tr[r 3 G(p,^ m )], (3.9) 

p P p,ui m 

where the renormalized single-particle Green's function G 

G(p,iu m ) = G (p,iuj m ) +G (p,iu m )E(p,iuj m )Go(p,iuJ m ), (3.10) 

involves the correction from the lowest order quasiparticle self-energy due to coupling with 
collective fluctuations 

E(p,iLO m ) = ^ ~T ~ — + UE u (q,iu n ))T + G (p + q,iw m + iw„)r_ 



+ (1 + £/H 22 (q, iu n ))T-G Q (p + q, iu m + iu n )r + 

- 2UE 12 (q,iv n )T + Go(p + q,iuj rn + iuj n )T + . (3.11) 



Here, ^(q, iv n ) = det[l + Z7H(q, iv n )], and t± = (ri ±zt 2 )/2. In obtaining this result, we 
have carried out the //-derivative on By in Eq. ()2.11|) by using the identity, 

^ = -G r 3 G . (3.12) 
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For example, 
djjL 



i [Tr[r 3 G (p, iw m )T-G{p + q, iu m + iu n )r + G (p, iu m )} 







Pi*Wn 



+ Tr[r3Go(p + q,^ m + w„)r + G(p,^ m )r_Go(p + q,^ m + w„)]]. (3.13) 
In the presence of a supercurrent, the Green's function analogous to Eq. ()3.10|) is given 

by 

g(p, iu m ) = g Q (p, iuj m ) + g (p, zcj m )E(p, iu m )g (p, ico m ), (3.14) 

where E is given by Eq. (j3.11j) but with g in Eq. ()3.7j) replacing G . Substituting Eq. 
(l3~T4l into Eq. we obtain 



Pn = Pn+ Pn- 



(3.15) 



The Fermi contribution is given by Eq. (J3.8|) . The bosonic fluctuation contribution p 1 ^ is 
given by 



9 B n 



2 d 

g ^^- Tr ^o(P,^m)E(p,ZU; m )^o(P,^m)]Q^O 



<9 



Z 3 P 7t ^Q, 
where we have used the identity 



Tr 



£(p,iu; m )- 



dgo 
diujm 



-9o9o- 



Substituting Eq. (|3. 1 lj) [with G — > g ] into Eq. ()3.16|) . and using the identity 

dg dg 



diujr, 



(3.16) 



(3.17) 



(3.18) 



we find 



B _2_ v _d_ v z U 



x Tr 



[1 + UE. 1 i{q 3 iv n )]T + g (jp + ,iu m + iu n )T-- 



(p-,iu m ) 



+ [1 + UE 22 (q, iv n )]r^g (p + , iu m + iu n )r+ 



diu ri 
dg (p„,iu n 



- 2C/5 12 (q, iv n )T + g (p + , iu m + iv n )r + 
2m d . , 



diu n 
dg (p-,iu m ) 



ditUr, 



[1 + Us. n (q,iu n )\- 



+ [1 + £7n 22 (q, «z/ n )J — 2U^ V2 (%iv r 



(3.19) 
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where p± = p ± q/2. We note that the correlation functions Hi 3 -(q, ii/ n ) appearing in Eq. 
(j3.19|) are defined as in Eq. (jA3|) in terms of the single-particle Green's functions go m the 
presence of a supercurrent. In Eq. (|3.19|) . the Q 2 -derivative only acts on Q z in the Doppler 
shift term a p = Q ■ p/2m in g in Eq. (J3.7j) . It does not act on the Q z in the shifted chemical 
potential p = p — Q 2 /8m. 

To summarize, we have shown that the total normal fluid density p n associated with 
fermionic and bosonic degrees of freedom is given by the sum of their contributions: 

n , tt- i ■ ^dZ u (q,iis n ) . dE 12 (q, iu n ) i 

+ [l + l/^Cq,^)]—^- 2^ 12 (q,^ w ) — j Q ^ Q . (3.20) 

We note that the superfluid density p s can be also obtained from a current correlation 
function The present derivation based on calculating the change in the single- 

particle Green's function from Eqs. (|3.2|) and (|3.14|) is equivalent to calculating the current 
correlation function to first order in v s , taking into account both self-energy and current 
vertex corrections. In Eq. (j3.2j) . the (^-derivative of the first term in Eq. (j3.14j) gives the 
bare current response function. The Q z -derivative of (?o in the second term in Eq. ()3.14|) 
gives the self-energy corrections, while the Q 2 -derivative of E gives the vertex corrections. 
See also Appendix A of Ref. [lOj for further discussion. 



IV. NUMERICAL RESULTS FOR SUPERFLUID DENSITY 

In this section, we present numerical results for the superfluid density p s , starting from 
the expression given in Eq. ()3.2()|) in Sec. III. As we have noted earlier, the expression for p s 
derived in Ref. [l^l would give identical results. We emphasize that these numerical results 
for p s use the renormalized values of both A and p{4] which determine the BCS quasiparticle 
spectrum over the entire BCS-BEC crossover. 

Figure |Hl shows the calculated superfluid density p s in the BCS-BEC crossover. The 
spurious first-order behavior near T c in the strong-coupling regime is also seen in the self- 
consistent solutions for A and p in Fig. El As discussed in Sec. II, this behavior near T c 
is removed through a more sophisticated treatment of fluctuations, which could lead to the 
correct second order phase transition. We plot our calculated results for p s close to T c , in 
spite of this problem. In particular, we note that the predicted value of T c in the NSR theory 
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FIG. 6: Calculated superfluid density p s in the BCS-BEC crossover. The self-consistent solutions 
for A and \l shown in Fig. [21 are used. The dashed line shows T c [see Fig. EI a)]. 

of the BCS-BEC crossover is a good approximation . The NSR bendover is much less in 
evidence in the case of a narrow Feshbach resonance, as considered in Ref. 0]. In Appendix 
C, we prove analytically that our expression in Eq. ()3.20|) gives p s = n at T = and also 
that p s vanishes as A — > (normal phase). Getting these two limits correctly (as shown in 
Fig. EJ) is very important in any theory of the superfluid density. 

Figure HI shows p s as a function of temperature in the BCS regime, unitarity limit, and 
the BEC regime. We note that p s in the BEC regime [(kpa s ) ~ x = 2] is in good agreement 
with the superfluid density p s of a weakly interacting gas of N/2 Bose molecules described 
by Bogoliubov-Popov excitations in Eq. ()2.17|) . as one expects in the extreme BEC limit. 
More precisely, in this limit one can show (see Ref. for details) that 



where is given by the Landau formula for the normal fluid of an interacting Bose gas 



Ps 



n- p n 



~ n — p. 



B 



(4.1) 



n ■ 



Pn 




2 dn B (uj (l ) 
duJd 



(4.2) 
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FIG. 7: Superfluid density p s as a function of temperature in the BCS region (solid circles), 
unitarity limit (solid triangles) and BEC regime (open circles). 'BCS' labels the mean-field BCS 
result, given by p s = n — p^ with p = ep- 'BEC gives p s for a dilute Bose gas with N/2 bosons 
described by the excitation spectrum in Eq. ()2.17j) . 

Here, M = 2m is the Cooper-pair mass and ns(c<; q ) is the Bose distribution function. The 
excitation energy u; q is given by Eq. ()2.17|) . calculated with the correct values of A and p,. 
The curve labeled by BEC in Fig. corresponds to the result obtained using Eqs. (j4.1|) 
and ()4.2|) . This shows the importance of a consistent treatment of fluctuation effects in 
calculating p s , A, and \x. 

As one approaches the weak-coupling BCS regime, pair fluctuations become weak, so that 
in this limit the fermionic contribution p 1 ^ becomes dominant. Thus one has 



where p^ is given by the Landau formula for the normal fluid in Eq. (|3.8jl with the BCS 
quasiparticle energies E p . The curve labeled by BCS in Fig. [7| corresponds to the result 
obtained using Eqs. (|4.3jl and (|3.8j) . 

V. CONDENSATE FRACTION IN THE BCS-BEC CROSSOVER 

In this section, we calculate the condensate fraction N c in the BCS-BEC crossover. 
The cerate .action in the superihnd phase is ,no 8 t convenient* dehnedB 
as the maximum eigenvalue of the two-particle density matrix, p2(r, r', r", r'") = 



Ps 



n- pn 




(4.3) 
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(^{(r)^{(r , )^ i (r")VT(r w )) ) where ^ ff (r) is a fermion field operator. The condensate frac- 
tion N c is given as the maximum eigenvalue, of order N. When only one eigenvalue is 
O(N), we find 

p 2 (r, r', r", r'") = N c ^ (v, r')0 o (r", r'"), (5.1) 

where terms of order 0(1) have been ignored. Here 0o( r > r is the (normalized) two-particle 
eigenfunction of p 2 with the eigenvalue N c . The off-diagonal long range order of a Fermi 
superfluidfsl is characterized as, for large separation of (r, r') and (r",r w ), 

p 2 (r,r',r",r'") = <^(r)^(r')> (^(r> T (r'")>. (5.2) 

Comparing Eqs. (|5.1jl and (|5.2jl . the condensate fraction N c is seen to be the normalization 
factor of the Cooper-pair wavefunction, $(r, r') = (^(r)^^')) (see, for example, Ref. Q|), 

N c = J drdr'|$(r,r')| 2 . (5.3) 

In physical terms, the maximum eigenvalue N c describes the occupancy of two-particle states. 
In a uniform Fermi superfluid, the BCS mean-field approximation gives 

*(r, r) = E(4T C -Pl) e ^ (r_r ' ) = E T^r tanh § E p e^'\ (5.4) 
p P 2h P 2 

In the strong-coupling BEC regime (where p <C —Bp), we can set tanh (3E p /2 = 1 in Eq. 

fl5.4|) . In this case, substituting Eq. (|5.4jl into Eq. (|5.Hjl . the mean-field expression for the 

condensate fraction (= A^o) in the BEC regime reduces to 

A 2 A 2 
p p p ''P 

In obtaining this expression, we have used the fact that |p| ^> A in the BEC regime jl3^. 
More generally, in terms of the single-particle Green's functions, one can write Eq. ()5.3|) 

as 

Nc = 7p ^ G 2 i(p,iu} m )G 12 (p,iu}' m ). (5.6) 

To calculate the strong-coupling effects on N c , we substitute Eq. ()3.1()|) into Eq. (|5.fi|) . Since 
this Green's function G = Go + GoEGo only includes first order self-energy corrections, we 
only retain the correction terms to iV c to 0(E), giving 

N = N c0 + 5N C . (5.7) 
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Here, the mean-field component N c q is the BCS Fermi quasiparticle contribution 



iV c0 = - £ Gf{p,iu m )G 12 {p,iu' m ) 




(5.8) 



The first order fluctuation contribution 5N C is given by 



1 



[^(P. ^m)Tr[r_G (p, zu4)£(p, iu' m )G (p, iu' m )\ 



ff 1 



+ Tr[r + G (p,^ m )S(p,zu; m )G'o(p,^m)]Go 2 (p,^) . 



(5.9) 



The correction term <5iV c in Eq. (|5.9|) is not important in the weak-coupling BCS regime, 
where fluctuation effects clearly can be ignored. Figure El shows that 5N C is also negligi- 
bly small in the strong-coupling regime. Thus N c is well approximated by the mean-field 
expression in Eq. ()5.8|) over the entire BCS-BEC crossover, at least in our NSR-type ap- 
proximation. We recall that the same pair fluctuations made a large contribution to p s as 
we went from the BCS to BEC region. The difference is that by definition, 5N C in Eq. ()5.9|) 
arises from self-energy corrections to the single-particle anomalous Green's function G\2- 
There is no distinct bosonic contribution, such as in the normal fluid density. Thus it is 
not unexpected that the fluctuations are a small correction to N c . 

We note, however, since fluctuations in the Cooper channel are taken into account in 
the equation of state in Eq. ()2.15|) . they modify the condensate fraction N c given by Eq. 
(|5.8|) . For example, let us consider the BEC regime at T = 0, where the gap equation gives 
p = —l/2ma 2 s . In this limit, the number equation reduces to 



gives the quantum depletion from the molecular condensate due to the effective interaction 
Um between Cooper pairs. Recently, the mean- field result in Eq. (15 .8j) has been used to 
study the condensate fraction N c in a superfluid Fermi gas at T = OjUj. In this case, the gap 
and number equation in the mean-field approximation for the BCS-BEC crossover reduce to 



N = N c0 + N d , 



(5.10) 



where N c0 is given by Eq. ()2.19|) and [taking the T — ► limit of Eq. ([2.16)1 ] 




(5.11) 




(5.12) 



19 



0.5 



(kfa s y J = 2 



0.1 



0.4 - 



0.3 



0.2 - 












0.05 



0.1 



0.15 



0.2 



0.25 



T/£ F 



FIG. 8: Condensate fraction N c in the strong-coupling BEC regime. The solid line shows iV c o and 
the dashed line includes the small correction 5N C from self-energies due to pairing fluctuations. 



In the BEC regime, while Eq. (|5.12jl again gives \i = —l/2ma%, Eq. (jo.lHj) reduces to 
N/2 = N c0 in Eq. (|2.16|) . As expected, the depletion N d at T = from the condensate 
due to the interaction between Cooper pairs is omitted when the BCS-BEC crossover is 
described by mean-field approximation . 

Figure 01 shows the condensate fraction in the BCS-BEC crossover regime at T = 0. 
Because of the omission of the quantum depletion effect, the simple mean-field result (MF) 
given in Ref. [ill is larger than our result (N cQ ). In the region (k-pag)' 1 > 1.5, our results 
are well described by the condensate fraction for a superfluid molecular Bose gas given by 
Eqs. (loTUl) and (I5~TT| (labeled BEC in Fig. 

Figure 03 also compares our T = results with those obtained by quantum Monte Carlo 
(MC) simulations [3]. The latter calculation gives results consistent with clm = 0.6a s j3|. In 
contrast, our NSR theory gives the larger mean-field molecular scattering length au = 2a s . 
As a result, we overestimate the magnitude of the depletion and thus our values for iV c are 
smaller than the MC calculation in the BCS-BEC crossover regime. The measurement of 
the depletion deep in the BEC regime would be a useful way of determining the magnitude 
of a M . 

Figure ITU1 shows the condensate fraction in the BCS-BEC crossover at finite temperatures. 
In the weak-coupling BCS regime, the condensate fraction N c is very small even far below 



p 




(5.13) 
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FIG. 9: Calculated condensate fraction N c in the BCS-BEC crossover at T = (solid line). In 
this figure, as well as in Figs. HUI and fTTl we only show N c q. 'MF' shows the condensate fraction in 
the case when A and fi are determined in the mean-field results in Eqs. (|5.12jl and (|5.13|) . 'BEC 
shows the result for a Bose gas described by Eqs. (|5. 1U|) and (|5.11|) . The solid circles shows recent 
Monte Carlo results for iV c |l2| for comparison. 

T c , because only atoms very close to the Fermi surface form Cooper pairs which are Bose- 
condensed. In this regime, Fig. ^2 shows that the temperature dependence of N c is very 
well described by the weak-coupling BCS result. In the crossover region, the temperature 
dependence of N c deviates from the simple BCS result, as shown by the case (/cpa s ) _1 = in 
Fig. ^2 in the BEC limit, Fig. ITU1 shows that the condensate fraction at T = approaches 
N/2, reflecting the fact that all atoms form Cooper pairs which are Bose-condensed. In this 
regime, Fig. ^2 shows that the temperature dependence of iV c agrees with the condensate 
fraction for a dilute Bose gas in the Popov approximation given by Eq. (|2.1fi|) . Figure El 
shows that the BEC picture is a very good approximation when (kpa s ) _1 > 1.5. The fact 
that N c agrees with the Popov theory for a weakly interacting molecular Bose gas shows 
that the superfluid phase transition in this regime is dominated by the thermal depletion of 
Cooper pair condensate, and not by the dissociation of Cooper pairs characteristic of the 
weak-coupling BCS regime. 
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FIG. 10: Condensate fraction N c as a function of temperature in the BCS-BEC crossover. In this 
calculation, the self-consistent solutions for A and p, in Fig. |2] are used. 

VI. CONCLUSIONS 

In this paper, we have calculated the superfluid density p s and condensate fraction N c 
in the BCS-BEC crossover regime of a uniform superfluid Fermi gas at finite temperatures. 
We have included strong-coupling fluctuation effects on both p s and N c within a Gaussian 
approximation. The same fluctuation effects were also taken into account in calculating the 
superfluid order parameter A and Fermi chemical potential p in the BCS-BEC crossover, 
within the NSR theory 3, H, Q - 

The expression we used to calculate p s was derived from the single-particle Green's func- 
tion in the presence of supercurrent as given by Eq. (j3.1|) . which brings in self-energy 
corrections due to dynamic pair fluctuations. In this paper, we have concentrated on the 
explicit numerical calculation of p s within a Gaussian approximation. In contrast, our com- 
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FIG. 11: Condensate fraction N c as a function of temperature in the BCS regime, unitarity limit, 
and the BEC regime. The curve 'BCS' is the weakly-coupling BCS result. The curve labeled 'BEC 
is the condensate fraction for a Bose superfluid determined by Eq. (|2.16j) . Results are normalized 
to values at T = 0. 
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panion paperjiyj uses a different (but equivalent) formulation which exhibits the structure 
of p s in a more direct fashion, in particular, the relation to collective modes. 

Our result for the normal fluid density in Eq. ()3.2U|) naturally separates into a mean- 
field part associated with fermions p^ and a bosonic pairing fluctuation contribution p^. As 
discussed in Ref. jlO | . is given by the Landau excitation formula in Eq. (|3.8|) in the whole 
BCS-BEC crossover. In the strong-coupling BEC regime, p^ is negligible and p% reduces to 
the Landau formula for the normal fluid of a Bose gas of tightly bound Copper pairs jlo|. 
However, in the region near unitarity, p 1 ^ is not expected to be given by a Landau-type 
formula because the bosonic pairing fluctuations are strongly damped. It is in this region 
that the numerical calculations for p s reported in this paper are especially useful. 

The superfluid density is a fundamental quantity in two- fluid hydro dynamics [22] and we 
will use our results in a future study of hydrodynamic modes in the BCS-BEC crossover 
regime of a Fermi superfluid at finite temperatures. 

In contrast to the superfluid density, the mean-field expression for the condensate fraction 
iV c is a good approximation even in the strong-coupling BEC regime. The fluctuation con- 
tribution to N c gives rise to the non-condensate component. In the BEC regime, we showed 
that the fluctuation contribution gives the condensate depletion iV d due to the effective 
interaction Um between Cooper pairs, which is finite even at T = 0. 



23 



In the BEC regime, the strong coupling theory presented in this paper reduces to that 
of a weakly interacting Bose gas of molecules, with an excitation spectrum given by the 
Bogoliubov-Popov approximation. This is also the origin of the spurious first-order phase 
transition our theory exhibits (see Figs. E]and[7J). This is a well-known problem in deal- 
ing with dilute Bose gases 25]. The recovery of the second order phase transition in the 
entire BCS-BEC crossover, as well as the normalized magnitude of the effective interaction 
between Cooper pairs 



, both require the inclusion of higher order fluctuations past the 
NSR Gaussian approximation which we have used. In this regard, we have emphasized that 
calculating the value of p s is very dependent on using the strong-coupling approximation for 
A and p as well, quantities which determine the single-particle excitation spectrum. Thus, 
when we calculate A and p beyond the Gaussian fluctuation level, we also need to improve 
the microscopic model used to calculate p s . The approach presented in this paper, as well 



as in Ref. 



10(| , can be the starting point for such improved calculations. 
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APPENDIX A: NUMBER EQUATION IN THE BCS-BEC CROSSOVER 

We briefly review some of the formalism developed in Ref. Q|. The thermodynamic 
potential Q = Q F + Q B consists of the mean-field part fl F due to Fermi quasiparticles and 
the correction term Q B due to boson fluctuations. The mean-field contribution Q F is given 
byflEil 



V F = f + - E p ] - \ £ ln[l + e"/"H (Al) 
17 P P P 

In the NSR Gaussian approximation, the fluctuation contribution is the sum of diagrams 

fl B shown in Fig. ^ namely 



1 



where 



n° 1 + n° 2 + i(n° 2 -iT 



^ rro 
V 



n° - n° 

LL U ii 22 



ii 2: 
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u° u + n° 2 - i(n? 2 - n° 21 ) j 



The correlation functions 11° ■ are given within the mean-field approximation by 



n°(q>n) = -^dre^(T T K q (r)a- q (0)}} 



/3 



£ T-iGolP + q/2, io; n + iv n )TjG (p - q/2, iw r , 



P,^n 



(A2) 



(A3) 



(A4) 



and II22 describe the amplitude fluctuations and phase fluctuations of the order param- 
eter, respectively, while II^ and 11^ are the coupling of these fluctuations. Doing the Fermi 
Matsubara frequency sum over ui m , one finds J, 

n?i(q,w n ) 



18] 



yv ( 1 _ ^p+q/2^p-q/2 ~ A 2 \ 
p \ -£p+q/2-Ep-q/2 / 



-^p+q/2 - -£p-q/2 

(£ p+q/2 - E p _ q/2 )2 + vl 



A 2 



p+q/2 + -frp-q/2 



x [f(E p+(l/2 ) - /(^p- q/2 )] 

_ Wl+ ^p+q/ 2 ^p-q/ 2 

p V #p+q/2-Ep-q/2 J (-^p+q/2 + -Ep-q/2) 2 + V% 

x [i-/(£; p+q/2 )-/(£; p _ q/2 )], 



(A5) 



i4(q>„) = Efi- ep ; q/2 ^ /2 + A 1 

p V -^p+q/2-^p-q/2 / 



-^p+q/2 ~ -^p-q/2 
(Ep+d/2 - ^p-q/2) 2 + Vl 
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X [/(£p +q /2) - /(^p-q/2)] 

£p+q/2£p-q/2 + A 2 \ -^p+q/2 + ^p-q/2 




-Ep+q/2-Ep-q/2 / (-^p+q/2 + -Ep-q/2) 2 + ^, 

"/(^p-q/2)], (A6) 



no / \ i ^p+q/ 2 ^p-q/2 



P 



V-^p+q/2 -^p-q/2/ (-^p+q/2 - -Ep-q/ 2 ) 2 + ^ 
X [/(^p+q/2) " /(^p-q/2)] 
/ £p+q/2 £p-q/2 \ 



^ V^P+q/2 ' Ep-l/z) (-^P+q/2 + -Ep-q/2) 2 + V\ 
X [1 - /(£ p+q/2 ) - /(£p- q/2 )] 

= -^(q,^). (A7) 

We note that for zero superfTuid flow (v s = 0), the matrix elements My of the inverse 
fluctuation propagator discussed in Ref. jlCJ are related to H^- by 

+ (A8) 

In contrast to Ref. jlo| . Eq. (jA3|) splits the fluctuations explicitly into phase and amplitude 
components. 

In calculating AT = —dfl/dfi, we note that, in principle, both £ p = e p — \i and A depend 
on [L. Thus, one needs to calculate dVL/dA = dVt F /dA + dVt B /OA. The first term d£l F /dA 
vanishes by definition, since the gap equation ()2.9|) is the saddle point solutions]. As 
noted in our previous paper j^], strictly speaking, the second term dfl B /OA is a higher order 
correction within the NSR Gaussian treatment of fluctuations. Thus in a consistent theory 
built on treating fluctuations to quadratic order, we should omit the term dfl B /dA. Within 
this approximation, we obtain the expression for A" given in Eq. (|2.11)1 . If we retained 
terms involving dQ B /OA in Eq. (|2.3|) . the resulting equation for A" would not reduce to the 
correct result at T c in the limit A — ► 0. On the other hand, the importance of this class o 



21 



33| . In particular, Ref. 
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higher order correction terms has been pointed out in Refs. 
has shown that the correct molecular scattering length in the BEC regime (cim = O.Qa s ) is 
obtained when this term is included. See Ref. il(J for additional discussion. 
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APPENDIX B: RELATION TO THE COUPLED FERMION-BOSON MODEL 



In this appendix, we compare the results for a single- channel model given by Eqs. ()2.14j) 
and (J2.15|) with our previous results for the coupled fermion-boson (CFB) model ^J. The 
CFB two-channel model involves treating the Feshbach resonance explicitly by the Hamil- 
tonian, 

#CFB = ~ ^H-Cpa ~ U hg C P+qT C P'-qi C P'l C Pt 

P' a p>p',q 

+ J2i4 + 2U - MKK + 3rE[^W,l C -p+q/2,t + ( B1 ) 

q p,q 
Here, 6 q is the creation operator of a molecular boson associated with the Feshbach res- 
onance, with the kinetic energy = q 2 /4M. The threshold energy 2v of the Feshbach 
resonance can be tuned by adjusting a small external magnetic field. g T describes the Fesh- 
bach coupling between atoms and a molecule. Uf, g is a nonresonant weak-interaction, which 
is taken to be attractive in Eq. (jBlj) . Since one molecule is a bound state of two Fermi 
atoms (in different hyperfine states), we take M = 2m and must impose a conservation law 
for the total number of Fermi atoms. The latter constraint has been taken into account in 
Eq. (IBljl by taking the chemical potential for the molecules to be 2\x. 

In a previous paper 4], we used this two-channel CFB model to treat the superfluid 
properties in the BCS-BEC crossover region within the NSR theory. The resulting coupled 
equations corresponding to Eqs. (j2.9J) and (j2.11|) are given by 



1 = E T^r tanh f E p , (B2) 
P 2-fcp 2 

N = 2<f> 2 m + N° F + 2N° B - i-A ]T In det \l + [U- <? r 2 £> Q (q, iu n )]E(q, iu n )} . (B3) 



2pdfi 



Here, U e s = U-^ g +g 2 j{2v— 2/i) is an effective pairing interaction associated with the Feshbach 
resonance. <p m = (6 q =o) is the BEC order parameter and § 2 m is the number of Bose-condensed 
Feshbach molecules. Because of the resonance coupling between atoms and molecules, <p m 

n 

is related to the BCS order parameter A by the relation <p m = —g r A/U(2u — 2/i)|4|. = 
Eq Wfi(^q + 2v — 2/i) gives the thermal occupation of the non-condensed molecules, where 
ns{x) is the Bose distribution function. -D°(q, iv n )~ l = i v n T z ~ (£q + 2z/ — 2/i) is the 2x2- 
matrix single-particle Green's function for a free Bose gas of molecules. In Eq. (jB2|) . the 
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BCS single-particle excitation spectrum E p is given by E p = yg* + A 2 , where £ p = e p — fi 
and A = A — g r (p m is the composite order parameter involving the A and (j> m coupled order 
parameters. For more details, see Ref. j^J. 

Let use now consider the case of a broad Feshbach resonance, where the coupling g r is 
very large (g T \/n ^> £f)- hi this limit, the effective pairing interaction U e s defined below Eq. 
(IB3J) can be strong even when the threshold energy 2v is still much larger than the chemical 
potential \i ( < £p). In this case, one can neglect /i in U e s in the interesting BCS-BEC 
crossover regime. Since 2v is the lowest excitation energy of Feshbach molecules, we can also 
neglect 20^ and 2N^ in Eq. ()B3|) when 2v ^> 2ep. Similarly, for 2u ^> 2ep, dynamical effects 
of Feshbach molecules are not important, and we can use the approximation -D°(q, iu n ) — 
D°(0, 0) = — l/(2u — 2/i) ~ —l/2u in Eq. ()B3|) . The end result is that the coupled equations 
(IB2jl and (|B3)1 for the two-channel model reduce to a form analogous to Eqs. ()2.9j) and ()2.11|) 
in the single-channel model. We need only replace A by A|3J] and U by the two-particle 
effective interaction associated with the Feshbach resonance, 

U% = U be + 9 f v (B4) 

Thus we see how the single-channel description of the BCS-BEC crossover described by 
Eqs. (|2.9J) and (j2.11|) is a special case of the two-channel model in the case of a broad 
Feshbach resonance 2v ^> 2ep. The single-channel scattering length a s is related to the 
two-channel CFB model parameters by 

^ = -K - 4. (B5) 
m s 2v y ' 

Here, U^, g , g r , and v are all renormalized quantities, given bvj^ 

TT V Wc — 

Uhe = ; — i , q r = 1— , 2v = 2v — q r , . Bo 

The renormalized unitarity limit (a s — ► ±oo) corresponds to 2z> = 0. 

APPENDIX C: ANALYTICAL RESULTS AT T = AND A -> 

We first prove that the normal fluid density as given by the expression in Eq. (13.20)1 
vanishes at T = 0. At T = 0, clearly the first term vanishes. For the fluctuation 
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part we find from explicit calculations that one can change the Q 2 -derivative into the 
^-derivative. Then, we find 



Pn(T = 0) 



1 + UE(q,iis n ) 



-- Y -2-—. rr In det 

2 q 2 f°° r d 2 

rE^ / dzn B (z)lm — In det [1 + UE(q, %v n — >• z + iS)] 



7T q 2m 



d 



-Y^-lm — lndet[l + [/H(q,z + zS) 
7r q 2m L az 



z=0 



(CI) 



In obtaining the last line, we have used that Ub(z) = —9{—z) at T = 0. One can easily show 
that Eq. ()C1|) vanishes, by noting that E^{z — > — oo) = and n n (z = 0), II 22 (z = 0), and 
= 0)n 2 i(2; = 0) are all real quantities. Thus, we have shown that in our microscopic 
model, p n = p^ + Pn — at T = 0, or p s = n. 

We next show that p s vanishes in the limit A — ► 0, i.e., in the normal phase above T c . To 
see this, we note that the mean-field part p^ reduces to the number of Fermi atoms when 
A — > 0, namely, 

pf(r c )= 2 (C2) 
p 

For the fluctuation contribution p%, since H 12 = S 2 i = gl 2 = g 2 } = in the limit A — > 0, 
Eq. (I3.19|) reduces to 

U 



Pn(T c 



--E— Ep* 

[1 + UE n (q,iu n )]g 22 (p + ,iuj m + 



+ [1 + £/H 22 (q, iv n )]g (p+, iuj m + w n 



diuj m 
dgi 2 (p-,iuJm) 



(C3) 



Carrying out the Qz-derivative, we obtain 
1 



Pn(T c ) 



E E ZT^TT f[l + UE 11 (q, iu n )]Gf(p + ,iu; m + iv n ) ^"(P-, ^, 



p,u™ q,i/ n VyQjWn 



[1 + t/H 22 (q, ii/„)]G (p + , iu m + w n 



diuOr. 



dlUJr, 



x 



G (p+, 2w m + «z/ n ) — G (p_, iu r , 



9Gj 1 (p + ,za; m + zi/ n 



(C4) 
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Noting that rj(q,iu) = [1 - UE n (q, iv n )][l - UE 2 2(q, iis n )} and E n (q,iv n ) = -U(q,iu n ) 
when A = 0, and using the identity, 

1 vMr'H/' • i • ,dGf(p-,iuj m ) 22 . 9G J J 1 (p+, iu m + iu n 
- EJG ( P+ ,- m + ^)^— G (p_,«, m ) — 

9 1 

= |-n(q> n ), (C5) 

one may reduce Eq. (|C4|) to 

(r e ) = -| £ Yu M 1 - tm(q, «v B )]. (C6) 



The sum of Eqs. (IC2|) and ()C6|) equals the total number of Fermi atoms n as given in Eq. 
(12. 3 j) . Therefore, we have proven that p n = n, or p s = in the limit A — > (normal phase). 
This is an important requirement of any theory of p s . 
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